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I critically review a methodology of using clusters of galaxies as cosmological probes. The 
understanding of the abundances and spatial correlations of dark matter halos has been sig- 
nificantly advanced especially for a last few years. Nevertheless such dark matter halos are 
not necessarily identical to clusters of galaxies in the universe. This is a quite obvious fact, 
but it seems that the resulting systematic errors, on the cosmological parameter estimates for 
instance, are often neglected in the literature. 

1 Why clusters ? 

This is probably an easy question to answer. In fact, there are several reasons why clusters 
of galaxies are regarded as useful probes of cosmology, including (i) since dynamical time-scale 
of clusters is comparable to the age of the universe, they should retain the cosmological initial 
condition fairly faithfully, (ii) clusters can be observed in various bands including optical, X-ray, 
radio, mm and submm bands, and in fact several on-going projects aim to make extensive surveys 
and detailed imaging/spectroscopic observations of clusters, (iii) to the first order, clusters are 
well approximated as a system of dark matter, gas and galaxies, and thus theoretically well- 
defined and relatively well-understood, at least compared with galaxies themselves, and (iv) on 
average one can observe a higher- z universe with clusters than with galaxies. It is established that 
X-ray observations are particularly suited for the study of clusters since the X-ray emissivity is 
proportional to n\ and thus less sensitive to the projection contamination which has been known 
to be a serious problem in their identifications with the optical data. Also the recent progress of 
interferometric mapping technique of clusters via the Sunyaev-Zel'dovich effect enables one to 
observe the high-redshift clusters without suffering from the cosmological dimming oc (1 + z)~ 4 . 

For those reasons, there are many, perhaps already too many, papers which discuss cos- 
mological implications of clusters of galaxies from various statistical methods. Table |l| lists 
representative topics of cosmological studies of clusters, with some possible bias to my personal 
interest, and thus in particular the reference list is far from complete. 

In what follows, I discuss cosmological implications of the Sunyaev - Zel'dovich effect (§Q) 
and cluster abundances (§|3|), paying attention to the limitations and possible systematics of 
current theoretical modeling of galaxy clusters. 



Table 1: An incomplete list of "cosmology with clusters" 



topic 


quantities 


references 


distance indicator 


Ho, Qq, Ao 


1 - 5 


peculiar velocity field 


^pec 


6- 9 


mass, temperature and luminosity functions 


08, flo 


10 - 20 


spatial clustering and its evolution 


£(r,z), b(r,z) 


21 - 29 


baryon fraction and dark matter 


Qb, ^0 


30 - 33 


CMB anisotropy through the SZ effect 


5T/T 


34 - 36 


universal density profile/nonlinear clustering 


p(r), P(k) 


37 - 45 



2 Distance and peculiar velocity from the Sunyaev-Zel'dovich effect 



There have been many attempts to determine the Hubble constant Hq (= lOO/i kms -1 Mi 
from the Sunyaev-Zel'dovich (SZ) temperature decrement and the X-ray measurements" 
fact, the observation of the SZ effect determines the angular diameter distance d\(z) to the 
redshift z of each clusten'B. For z <C 1, d\(z) is basically given only in terms of Hq (~ cz/Hq, 
with c being the light velocity). If z ^ 0.1, however, the density parameter Qq, the dimensionless 
cosmological constant Ao, and possibly the degree of the inhomogeneities in the light path make 
significant contribution to the value of dx(z) as well. As far as I understand, Figure p] is the first 
Hubble diagram from the SZ effect published in the literatureP Note that the data for CL0016+16 
plotted in this figure have been corrected later. Figure |l| clearly indicates the potential difficulty 
of determining the cosmological parameters from the SZ effect of an individual cluster. Of 
course this is not surprising because the estimate of the angular diameter distance from the 
SZ effect crucially rests on several idealistic assumptions including the spherical symmetry, 
isothermal temperature distribution, and no dumpiness of the gas. It is quite ironical that 
one can immediately point out the possible systematics by abandoning one of the assumptions 
because the SZ effect is based on the well-defined physics rather than on the unjustifiable, and 
thus unfalsefiable, empirical result which is the case for the other distance indicators. 




Figure 1: The angular diameter distance c£a estimated from the SZ effect as a function of z for seven clusters 
(symbols). Lines indicate theoretical curves for several sets of cosmological parameters (h, f2o, Ao). From reference 

4. 

Therefore it is essential to accumulate a statistically meaningful number of the SZ clusters. 
For this purpose in mind, we performed Smoothed Particle Hydrodynamical (SPH) simulations 
of clusters of galaxies in Cold Dark Matter (CDM) modelsB Figure || presents an example of a 
simulated cluster observed in different bands. This illustrates how an individual cluster exhibits 
the departure from the above idealistic assumptions. Figure [3| summarizes the distribution of 
the estimated Hubble constant from the projected X-ray surface brightness profile (left) and the 
3D density profile (middle) and of the estimated peculiar velocity of clusters (right)® The result 
suggests that the intrinsic scatter reflecting the different cluster gas state is fairly large, especially 
at high redshifts, and that at least a few tens of clusters should be observed to determine 
the parameters within 10 percent accuracy. Of course this conclusion is heavily dependent on 
the extent to which the simulated clusters faithfully reproduce the statistical properties of the 



observed cluster samples in the universe. To properly answer this question, we should wait for 
an extensive survey of clusters with high-angular resolution. Incidentally, our recent SZ image 
(cfwhm = 13" at 150 GHz) of the most luminous X-ray cluster RX J1347-1145 [z = 0.45) 
indeed revealed complex morphological structures of the cluster region, and therefore the 
departure from the the spherical isothermal /3-model for the clusters may be more significant 
than currently thought. 




Figure 2: Projected views of a simulated clusters at z = 1 and z ft! 0. A box of (4Mpc) 3 (in physical lengths) 
located at the center of each cluster is extracted. The X-ray emission-weighted temperature (Tx), X-ray surface 
brightness (Sx), and the SZ surface brightness at mm and submm bands (|A7 mm | and AJ su bmm) are plotted on 
the projected X-Y plane by integrating over the line-of-sight direction (Z). The X and Y coordinates are in the 
physical lengths at the corresponding redshift, and related to the angular coordinate 9 from the cluster center as 
6<1a(z). At z « 0, <Ia(z) can be replaced by the real distance to the cluster from the observer. From reference 9. 



3 Cluster abundance 



Cosmological implications of cluster ab unda nce are discusse 
ing X-ray temperature functionj^Stllllo ma ss functio nr 4 ^ 

lc 



ising a variety of statistics includ- 
23 velocity functioiEJ'LJ and X-ray 
luminosity function and log N - log S re] ation I 1 ilMM^H HlM In what follows, we specifically 
consider the theoretical modelliZHISEJ for X-ray log ./V - log S relation in order to illustrate the 
systematic effect in cosmological conclusions derived from cluster abundance. 

The number of clusters observed per unit solid angle with X-ray flux greater than S is 
predicted from 



N(> S) 



dz d 2 A (z) c 



dt 

dz 



/ dS {l + z) 3 n M {M,z] 
Js 



dM <iT gas dL band 



dT Bail dL h 



dS 



(1) 



where c is the speed of light, t is the cosmic time, d A is the angular diameter distance, T gas and 
^band ar e respectively the gas temperature and the band-limited absolute luminosity of clusters, 
and hm(M, z)dM is the comoving number density of virialized clusters of mass M ~ M + dM 
at redshift z (we use the Press-Schechter mass function). 
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(a) z=1.0 f HAD =0.85±0.20 (c) z=1.0 f HJD =0.73±0.17 



(e) z=1.0 f v =0.87±0.10 
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Figure 3: Distribution of /h,2d {left), /h,3D {middle), and {right) for all of the simulated clusters (nine in total) 
at z = 0.05 and 2 = 1.0 viewed from three different line-of-sight directions. Different patterns of the histogram 
correspond to different clusters. The mean and lcr statistical errors are quoted in each panel. From reference 9. 



Given the observed flux S in an X-ray energy band \E a ,E b \, the source luminosity L band at 
z in the corresponding band [E a (l + z),Ej,{1 + z)\ is written as 

L hand [E a (l + z),E b (l + z)]=4irdl(z)S[E a ,E b ], (2) 

where d h = (1 + z) 2 d A is the luminosity distance. We adopt the observed L bol — T gas relation 
parameterized by 

L ^ = L ^(-^)° \ 1 + z ) < lO^^ergsec" 1 . (3) 

We take L44 = 2.9, a = 3.4 and £ = as a fiducial set of parameters. Then we translate L bol (T sas ) 
into the band- limited luminosity L band [T gas , E\, E%] taking account of the thermal bremsstrahlung 
and the metal line emissions. 

Assuming that the intracluster gas is isothermal, its temperature T gaa is related to the total 
mass M by 

\xm v GM 
7 3r vir (M, z f Y 

5 ^ i+ ^(^) 1/3 (ist^) 2/3 «; /3 ^ w 

where \i is the mean molecular weight (we adopt \x = 0.59), and 7 is a fudge factor of order 
unity. The virial radius r v i r (M, Zf) of a cluster of mass M virialized at Zf is computed from 
A v ; r , the ratio of the mean cluster density to the mean density of the universe at that epoch. 
Further we assume that the gas temperature evolves as 

T gas (M,z f ,z) = (i±^yr gas (M,z f ), (5) 
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Figure 4: The SZ maps of RX J1347-1145. (a) the 21 GHz map: 6' x 6' (2.4 ft" 1 Mpc x 2.4 /i" 1 Mpc). (b) 
the 21 GHz map after subtracting the central point source, (c) the 150 GHz map: 1'9 x 1'9 (0.75 h~l Mpc x 
0.75 Mpc). (d) the 150 GHz map after subtracting the central point source (assuming the flux of 3.8 mjy) 

overlaid with the X-ray contours. From reference 48. 

taking account of the quiescent accretion of matter after the formation epoch 

The left panel of Figure || shows the predictions for the X-ray LogA-LogS relations in 
various CDM models. As is well known now, many models with appropriate sets of cosmological 
parameters account for the observed cluster abundance. Several examples of such models are 
listed in Table || which exhibit the almost indistinguishable predictions for the X-ray LogA- 
LogS" relations in the ROSAT band. 

The degeneracy among those viable cosmological models can be broken by observing wider 
(i.e., increasing the statistics), deeper (at higher redshifts) and/or in multi-wavelength bands. 
The middle and right panels in Figure |5| plot the cluster LogA-LogS" predictions in hard X-ray 
(due to the thermal bremsstrahlung) , and in submm (due to the SZ effect). Figure || illustrates 
the extent to which one can break the degeneracy between as and m CDM models (n = 1, 
h = 0.7) using the multi-band observational data. 

Similarly the redshift-distribution of cluster abundances can be a very powerful discriminator 
of the different cosmological models. Figure |j] exhibits the redshift evolution of the number of 
clusters in different bands. As expected, the evolutionary behavior strongly depends on the 
values of S7o an d erg; the fraction of low redshift clusters becomes larger for greater Qq and 
smaller ag. This indicates that one may be able to distinguish among these models merely by 



Table 2: CDM models consistent with the ROSAT X-ray Log N - Log S relation. 



Model 


^0 


Ao 


h 


0"8 


a 


7 


L03 


0.3 


0.7 


0.7 


1.04 


3.4 


1.2 


O045 


0.45 





0.7 


0.83 


3.4 


1.2 


El 


1.0 





0.5 


0.56 


3.4 


1.2 


L037 


0.3 


0.7 


0.7 


0.90 


3.4 


1.5 


LOIq 


0.1 


0.9 


0.7 


1.47 


2.7 


1.2 




Figure 5: The LogiV-LogS' relations of galaxy clusters for CDM models in (a) the soft X-ray (0.5-2.0 keV) band, 
(b) the hard X-ray (2-10 keV) band, and (c) the submm (0.85 mm) band. Lines represent the models listed in 

Table 2. From reference 20. 

determining the redshifts of clusters up to z ~ 0.2. 

Since the above agreement between model predictions and available observations is so re- 
markable, the resulting conclusions on the values of Oo and a% look robust. Figure ^ explores 
the possible systematic effects on fio and as by changing the model parameters which describe 
the cluster gas. For a realistic range of those parameters, the constraints change around (10- 
20) percent level. As a matter of fact, this statement is somewhat misleading. The crucial 
assumption underlying the cluster abundance modeling, I believe, is the one-to-one correspon- 
dence between the observed galaxy clusters and the dark matter halos. The former is defined via 
the optical luminosity of member galaxies and/or X-ray luminosity of the cluster gas, while the 
latter is defined basically from the spherical collapse model. More specifically, we have a variety 
of definitions for clusters; optically selected clusters (or the Abell clusters), X-ray selected clus- 
ters, SZ selected clusters, dark matter halos defined through the spherical collapse, and halos 
directly identified from large cosmological simulations using a variety of selection criteria. They 
should not be identical. Of course I agree that assuming the one-to-one correspondence among 
those species is a good approximation. My point here, however, is that the assumption may 
easily affect the derived values of the cosmological parameters more than the systematic effect 
presented in the above totally under this assumption. 

I do not think that this has not been considered seriously simply because the agreement 
between model predictions and available observations is satisfactory. Since current viable cos- 
mological models are specified by a set of many adjustable parameters (see example Table |2|), 
the agreement does not necessarily justify the underlying assumption. Thus it is dangerous to 
stop doubting the unlikely assumption because of the (apparent ?) success. In this respect, I am 
always impressed by the marked contrast with the case of the SZ effect as a distance indicator 
(§|2|), where the simple model predictions and the observations do not agree perfectly (with many 
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Figure 6: Contour maps on the fio-o"8 plane in (a) spatially flat (Ao = 1 — fio) CDM models, (b) open (Ao = 0) 
CDM models, (c) spatially flat CDM-like models with the fixed shape parameter (F = 0.25), and (d) open CDM- 
like models with T = 0.25. In all cases, h = 0.7, a = 3.4, and 7 = 1.2 are assumed. Shaded regions represent the 
la significance contours derived in KS97 from the soft X-ray (0.5-2 keV) FogiV-FogS' . Dotted and solid lines 
indicate the contours of the number of clusters greater than S per steradian (10 2 , 10 3 , 10 4 from bottom to top) 
with 5* = 10 -13 erg cm -2 s _1 in the hard X-ray (2-10 keV) band and with 5V = 50mJy in the submm (0.85 mm) 
band, respectively. Fhick dashed lines represent the COBE 4 year result. From reference 20. 

scatters) and everybody talks of the systematics carefully enough. 
4 Summary and conclusions 

Here I presented a brief review of cosmological implications of galaxy clusters, specifically con- 
sidering the Sunyaev-Zel'dovich effect and the cluster abundance in some detail. I believe that 
the clusters are quite important and useful probes of cosmology and in fact they already proved 
to be successful in many respects. If one would like to go further and to extract more strin- 
gent constraints on the cosmological models, however, the one-to-one correspondence between 
virialized halos and observed clusters, whatever they mean, should be critically examined. This 
assumption is a reasonable working hypothesis, but we need more quantitative justification or 
modification in order to improve the cosmology with clusters. 
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